Hydro dynamic approximation for the nonlinear response of a metal surface 
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We present semi-classical and quantized hydrodynamic models to obtain the quadratic electronic 
response of a plane-bounded electron gas. Explicit expressions for the dynamic image potential 
experienced by charged particles moving near a jellium surface are derived, up to third order in the 
projectile charge. These expressions are employed to compute the image potential at all distances 
outside the surface. Though nonlinear corrections are found to be more important far inside the 
solid than outside, our results indicate that the nonlinear image potential is enhanced with respect 
to the linear image potential by a factor that is for Al as large as ~ l.fS near the surface in the case 
of a stationary particle {v 0) with positive unit charge e. 

7f.45.Gm, 73.20.Mf, 79.20.Rf 

I. INTRODUCTION 

The electronic response of a metal surface to an exter- 
nal perturbation enters a great variety of important prob- 
lems in surface science.u Surface spectroscopies employ- 
ing electrons, photons, atoms, or ions all involve some 
kind of electronic excitation at the boundary of the sur- 
face. In particular, the interaction of charged particles 
with solids haSrjepresented an active field of basic and 
applied physicsp and recently a great amount of research 
has been focused on the case of slow {v < 0.5 vo, vq being 
the Bohr velocity) highly charged ions {Zi >> Ar\Zi e 
being the ion charge) moving near a solid surface . 13 For 
these projectiles, the parameter ZiVq/v is not small and 
first-order perturbation or, equivalentlv. linear-response 
theories are not, a priori, applicable.El In the case of 
charged particles moving inside a solid, nonlinear ef- 
fects have been found to be cpjiieial in the interpretation 
of energy-loss measurements.ErEJ Nevertheless, the elec- 
tronic response of metal surfaces to the presence of ex- 
ternal charged particles, which differs significantly from 
that in purely two- or three-dimensional systems,— had 
been described so far within linear response theory.E3 

A central quantity in the interpretation of ion-surface 
collisions is the so-called image potential, which repre- 
sents the interaction between the incoming charge and 
the polarization charge that it induces on the surface. 
In the case of a particle of charge Zie located at a dis- 
tance Zq far from the surface, into the vacuum, this poten- 
tial approaches the long-ranEe. classical Coulomb image 
potential V™^ = —Z^e^ /4:Zq£3 For smaller values of zq 
the image potential differs significantly from its classical 
limit, the deviation from the classical result increasing as 
Zq decreases JIj 

The classical image potential acting between a point 
classical charge and a metal surface may be regarded 
as or.iginja.ting in the coupling with the surface plasmon 
field.ll30 Dynamical corrections to the classical image 
potential have been discussed in the framework of linear 



response theory ,E3E2l and recoil effects have been treated 
by including the exchange of virtual epccitations between 
the external charge and the medium E2l Preliminary re- 
sults for nonlinear corrections to the image potential as- 
sociated to the quadratic response-iof solid surfaces have 
been reported only very recently.E3 

Theoretical approaches commonly used to describe the 
electronic response of jellium jSiirfaces can be classified 
as being either hydrodynamicnJ in nature or based on 
the so-called random-phase approximation (RPA).E3 Hy- 
drodynamic approaches are appealing because of their 
relative mathematical and computational simplicity, and 
have been used with great success in the, description of 
collective phenomena at metal surfaces£3 Within a hy- 
drodynamic model, one assumes that the collective mo- 
tion of the electron gas may be described in terms of the 
displacement of the electrons from their original uniform 
state, and the electron system is characterized by the 
electron density and a velocity field. These quantities 
are then obtained by solving the well-known nonlinear 
Bloch hydrodynamic equations. If one docs not include 
quantum properties of the electron system we refer to the 
hydrodynamic model as classical, and as semi-classical if 
quantum properties are introduced through the defini- 
tion of the internal energy density. If one quantizes the 
hydrodynamic riiamiltonian on the basis of the existing 
normal modes,Ej we have the so-called quantized hydro- 
dynamic model which allows to apply standard methods 
of many-body perturbation theory. 

In a previous workJSJ we used the quantized hydro- 
dynamic model to describe the quadratic response of a 
homogeneous electron gas. We derived expressions for 
the quadratic polarization induced by a moving charged 
particle, and demonstrated that thjcy coincide with the 
plasmon-pole (PP) approxinmtionES to the more accurate 
quadratic RPA polarizationEj 

In this paper, we first develop semi-classical and quan- 
tized hydrodynamic models to derive the quadratic elec- 
tronic response of a plane-bounded electron gas, and we 
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then focus on the evaluation of the nonhnear dynamic 
image potential experienced by charged particles moving 
parallel to a jellium surface. In sections II and III semi- 
classical and quantized nonlinear hydrodynamic models 
are presented, respectively. In section IV numerical cal- 
culations of both linear and quadratic contributions to 
the image potential are reported, as a function of the 
distance from the surface. In section V our conclusions 
are presented. 



II. SEMI-CLASSICAL HYDRODYNAMIC MODEL 

Take an inhomogeneous electron system embedded in 
a nciitralizing ionic background. In the hydrodynamic 
liipi|t,E2l the total energy of the system can be expressed 
a^ (we use atomic units throughout, i. e., h — rrie — 
e-1): 



H = - I drne(r,t)|VV(r,i)|2 



drn(r,t) V;xt(r,0 



i J drn{r,t)V{r,t) 
drGK(r,i)], (1) 



where irrotational flow has been assumed, i. e., u(r,t) = 
— V-0(r, i), u{r,t) being a velocity field, and retardation 
effects have been neglected. ne(r, i) = no(r) -I- n{r,t) 
is the total electron density, with n{r,t) representing 
the deviation from the equilibrium static density no(r). 
V{r,t) is the induced electric potential, 14xt(i", rep- 
resents the external perturbation, and G[ne{r,t)] rep- 
resents the exchange, correlation and internal kinetic 
energies of the electron system. We neglect exchange- 
correlation contributions to G [ne{r, i)],p4Mhich we approx- 
imate by the Thomas-Fermi functionaEJ 

GK(r,t)] = A(3,r2)2/3 ln,(r,t)f' . (2) 

From Eq. (|l|), the basic hydrodynamic equations can 
be derived, i. e., the continuity equation, 



and the Bernouilli equation. 
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r]ip + 11; + fi = - |V-0r - U + 



SG[r 



SUe 



(3) 



(4) 



which conserves both momentum and energy. Here, fi 
is a Lagrangian multiplier (a constant), and is a posi- 
tive number representing the internal friction of the elec- 
tron gas which would appear as a consequence of the 
interaction with excitations not included in this descrip- 
tion. If the external perturbation is generated by a 
charge density ncxt(r, t), then the total electric potential 
U{r, t) = V{r, t) + Vcxt{Tc, t) is obtained from the Poisson 
equation 



These are nonlinear equations, difficult to solve. 
Within perturbation theory, we expand the induced elec- 
tron density n(r, t) and the velocity potential ■(/'(r, t) in 
powers of the external perturbation. 



and 



n = ni+ 71-2 + 



(6) 



(7) 



respectively, assuming that uq >> ni » n2 » ... and 
-01 >> i}}2 » .... Eqs. (H) and (jj) are then expanded in 
powers of the external perturbation, and partial differen- 
tial equations for the various orders of n(r, t) and t) 
are derived. 



A. Linear Approximation 

Up to first order in the external perturbation one finds, 
after introduction of Eq. ^ into Eq. (jj), the linearized 
hydrodynamic equations 



and 



77 V'l + V'l = H ni 

no 



ni = Tio V ipi, 



where Ui{r,t) is obtained from 
V^C/i = -47r Kxt 



(8) 



(9) 



(10) 



and where /3 = y/T/Sqp, Qf ~ (STr^ng)^/^ being the 
Fermi momentum. Though this is the value of the hydro- 
dynamic speed /3 predicted with use of the Thomas- Fermi 
functional G [ne(r, t)\ of Eq. (^, the value /3 = ^fzjhqp 
is expected to be more appropriate when high frequencies 
of the order of the plasma frequency are involvedEjO 

We consider a classical charged particle moving with 
velocity v outside of a semi-infinite metallic medium, 
along a trajectory that is parallel to the surface, thereby 
approximately simulating the experimental conditions 
when the projectile approaches the surface at grazing in- 
cidence. Hence, we take the external charge density at 
r = (rii , z) to be given by the following expression: 



t(r,t) = Zi5(r|| - vi)(5(z-.2o), 



(11) 



the vacuum occupying the half-space 2; > 0, and zq be- 
ing the distance of the trajectory from the metal. After 
Fourier analyzing both in time and in ry , one finds the 
following linearized hydrodynamic equations with vari- 
ables (2;;q, cj), q being the wave vector parallel to the 
surface: 



V^JJ = -47r [next - n\ . 



(5) 



13^ 

-icj + 77) tjji = —Ui H ni 

no 



(12) 
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and 



- icjni = HQ (V'l - q'^ipi), 



(13) 



where the prime denotes the derivative with respect to z, 
[/i(z;q, w) is obtained from 



(14) 



and 



'^oxt(z; (i,uj) = 2tt Zi 5{uj - q • v) (5(2 - zq). (15) 

Assmxiing a sharp density profile at the surface, i. e., 
7io(z) — ?ioO(— z) (6(z) is the Heaviside step function), 
the appropriate solutions of Eqs. (|l2|) and ( p^ ) are: 

ni(z; q, c^) = no (A^ ^ - g^) e^-- ^ e(-z), (16) 

V'i(z;q,c^) = -ic^ (Aie'^<'-^+Bie«-) e(-z), (17) 

C/<(z;q,w) =cJ2^le'^■'-^+cJ(cJ + i7/)Sle■^^ (18) 
and 

U> (z; q, cj) = Ci e-« ^ + Zi — (5(cj - q • v) e"' l^"^"' . 

(19) 

Here, J7;j^'^ (z; q, cj) is the Fourier transformed electric 
potential for z (less than, greater than) zero, Wp = 
(47rno)^/^ is the so-called plasma frequency, and 



(20) 



The constants Ai, and Ci are evaluated from the 
boundary conditions at the surface: 



(21) 



i?, = ^., , (22) 



g [2cj(cJ + i77)-Cj2(l+g/Aq,^)] 



and 



^ ^ 4^2^! a^g(l-g/Aq,^)(5(^-q-v)e-^-° 
' 9 [2c.;(c^ + i77) -c^2(i + ^/Aq_^)] ' 

The induced electric potential Vi (r, t) is the difference 
between the total and external potentials Ui (r, t) and 
14xt (r, i) , respectively. The image potential is defined 
as half of the induced potential at the position of the 
projectile times the projectile charge, and one finds 



2 !< 



1/g - 1/Aq,q.v 

[2 q • V (q • V + iry) - a;2 (1+ g/Aq,q.v) 



(24) 



This agrees with the-result obtained using either the spec- 
ular refljExion (SR)Efl or the semi-classical infinite barrier 
(SCIB)Ej model of the surface, as long as the hydrody- 
namic dielectric response function of the bulk material 
is used in these models. In particular, for a stationary 
charged particle (v = 0), one finds 



= -Zl I dq '^ N,0 2g.o_ 
1 2 %-t-Aq,o 



(25) 



This agrees with the early result of Eguiluz. 

Within the hydrodynamic model, one can define two 
characteristic screening lengths A. For a stationary 
charged particle (v = 0), A = f^/i^p, and in the case 
of a swift charged particle {v » qp), A = v/tUp. In the 
limit zq oo, the distance zq of the projectile from the 
surface being much larger than the screening length, one 
finds the classical image potential 

Vr = -^. (26) 
4zo 



B. Quadratic Approximation 

Up to N-th order (N > 2) in the external perturbation, 
one finds the hydrodynamic equations: 

1 



m=l 

k=l 



nti(5-3s) 



k^2/3 
^0 



3k 



E 



ki ,k2,...,k]sj = 



EN v^N 




(27) 



and 



N-l 



where Vn (r, t) is obtained from 
V^Vn = 47r77,r 



(28) 



(29) 



After Fourier transforming in time and in r|| , we obtain 
for N = 2 the quadratic approximation 
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(z;q,uj) 



no 



■ n2 



6 n. 



2 "i 



and 



itj ^2 = ?^o('02 ~ '?^V'2) + [n-iV^Vi + Vni • V?/'i 





(30) 



(z;q,aj) 

(31) 



where the prime denotes the derivative with respect to 
z, and the Fourier transformed potential Yi{z\(\^uj) is 
obtained from 



(32) 



In particular, for a stationary (v — 0) charged particle 
Eqs. (pO) and ^ yield 



2 

^2 ^2 - 7-2 «!, 



(33) 



where V2(2;; q, lo) is still obtained from Eq. (|32|). Solving 
these equations with the assumption of a sharp density 
profile at the surface, quadratic contributions to the in- 
duced density and potential are found to be given by the 
following expressions: 

n2(z;q,c^)= e(-z) [Aae'^-"" 



enoJ (27r)2 



^(Ak.o+Ak-q.o) ■' 



(Ak,o + Ak-q,o)' - , . 

X ^ ^ Jk /k-q 



(Ak,o + Ak-q,o) - A 



^2<(^;q,^) 



4-77 



q,0 



^2 e^-'-^ + Bz 



(34) 



^ .47r /• d^k /k/k-qe(^''"+^''--"^^ 

+ ^1 c 



6no ] (27r)2 (a^^q + Ak-q,, 



and 



where 



F2>(z;q,c^) = C2e' 



/q = (Aq^O-'?) e-«^ 
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(35) 



(36) 



(37) 



and Aq.ij is obtained from Eq. (^0|). The constants A^,-, 
B2 and C2 are evaluated from the boundary conditions 
at the surface. 

Hence, we find the quadratic contribution to the image 
potential of a stationary (v — 0) charged particle to be 
given by the following expression: 



2tt 



dO 



/q /qi /q-qi 



(Aq,0 + Aq,,o+Aq_q,,o)' 



(38) 



9 being the angle between q and qi. In the case of a 
non-dispersive electron gas {(3 — Q), this contribution to 
the image potential vanishes. On the other hand, in the 
limit as zo ^ 00 only the low-momentum form of the 
integrand of Eq. (|3^) contributes to the integration, and 
we find in this limit 



dqi qi e 



ISttw^Jo ^0 

27r 



Numerical integration yields 

^3 0-82 

~ ^1 18 ^6.2^4' 

and, therefore (see Eq. (p^)), 

1.93 X 10-2 



'(q+qi)zo 



4 20 



1 + Zi 



[zo/rs 



O {Zl) 



(39) 



(40) 



(41) 



This is a reasonable approximation for the image poten- 
tial of a stationary charge, as long as the distance from 
the surface Zq is larger than the characteristic screening 



length /S/wp, 



1. e., for zq^^Ts oo, ao = ?i /(mge^) being 
the Bohr radius. We note that in the high-density limit 
(r^ 0) there is no quadratic contribution to the image 
potential of a stationary charge, while at metallic densi- 
ties (rs ~ 2 — 6) quadratic corrections might give rise to 
an image potential that is at zq y/r^ ao larger than the 
linear image potential by a factor as large as 1.05 — 1.3 in 
the case of a stationary particle with positive unit charge 
e {Zl = 1). 



III. QUANTIZED HYDRODYNAMIC MODEL 

Within a quantized hydrodynamic model of the elec- 
tron gas, we first expand the hamiltonian of Eq. (|l|) in 
powers of the induced electron density n{r,t). After in- 
troduction of the Thomas-Fermi functional of Eq. (|^) 
into Eq. dl^), up to third order one finds 



where 



TT 3 . 2\2/3 5/3 

^ " 10 ^ ' "° ' 



Hn= dr 



ino|V^P + /L^2_ 1 ^ 
z 2 no 2 



(42) 



(43) 



(44) 



= -Zf 



3 



6 7rw4 j„ 



dqq dqi qi 



Hi 



dr 



■nlW'l^ 



18 ng 



(45) 
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and 



where 



dvnVcy 



(46) 



Hq is the Thomas-Fermi ground state of the static un- 
perturbed electron system, Hq represents the hnear de- 
viation from the ground state, Hi appears as a conse- 
quence of the nonhnearity of the electron system, and 
i?ext represents the contribution to the hamiltonian from 
the coupling with the external charged particle. 

We consider, as in the previous section, a semi-infinite 
electron system embedded in a neutralizing ionic back- 
ground, assuming a sharp electron-density profile at the 
surface. For each value of q (the wave vector parallel 
to the surface) there exist both bulk and surface normal 
modes of oscillation with frequencies given by the follow- 
ing dispersion relations: 



and 



,1/2 



(47) 



(48) 



respectively, where /? represents the speed of propagation 
of hydrodynamic disturbances in the electron system. As 
in the previous section, we choose (3 = ^Ji/bqp- 

Now we follow Ref. ^ to quantize the hamiltonian of 
Eq. ( ^^ on the basis of the normal modes corresponding 
to Eqs. ( ^7|) and (|4^), which we shall refer after quan- 
tization as bulk and surface plasmons, respectively. We 
find 



Ho — H^ + (f , 



(49) 



H^ and Hq being free bulk and surface plasmon hamil- 
tonians, respectively: 

< = ^ E [1/2 + <p] «q,p(i) (50) 



q,p>0 



and 



(51) 



Here and A represent the normalization volume and 
the normalization area of the surface, respectively, and 
Qq,p(i) and bq(t) are Bose-Einstein operators that an- 
nihilate bulk and surface plasmons with wave vectors 
(q, p) and q, respectively. The quantized Hi hamilto- 
nian, which contains the quadratic electronic response 
of the electron system, will be consider below. For the 
hamiltonian containing the coupling between the exter- 
nal particle and either bulk or surface plasmon fields, one 
finds 



Hc^t — H" 



Hi 



(52) 



H 



B/S _ 



j dvp,^t{v,t)4>''/^{v,t), (53) 



(f)^/^(^r,t) representing operators corresponding to the 
scalar electric potential due to bulk/surface plasmons. 
Outside the metal (z > 0), 



q,p>0 



and 



(54) 



(55) 



Xq,p(i) and Xq(i) representing operators associated to 
the electron density induced by the excitation of bulk 
and surface plasmon fields, respectively, 



and 



Xq,pW = aij,p(i) + a_q,p(i) 



XqW =&UO+&-qW, 



(56) 



(57) 



and fqp(z) and f^{z) being bulk and surface coupling 
hmctions. 



[p4+p2(^2 +^2/^2) +^4/(4^4)] 



1/2 



(58) 



and 



[q (g + 27q)]'/' 



(59) 



Here 7q represents the so-called invp*se decay length of 
surface plasmon charge fluctuations,E3 



_ 1 



-I3q- 



f3^q^ 



(60) 



In the absence of electron- gas dispersion (/? 0), the 
scalar electric potential ip-^{r,t) due to bulk plasmons 
vanishes outside the surface; hence, in this case probes 
exterior to the solid can only generate surface excitations. 

We derive the potential induced by the presence of the 
external perturbing charge as the|-expectation value of 
the total scalar potential operatorj^ll 



Vir,t) = 



<*o|*o> 



(61) 



where I^^/q > is the Heisenberg ground state of the inter- 
acting system, and where (p^{r,t) and (/)^(r,t) are the 
operators of Eqs. js^ ) and ( p5| ) in the Heisenberg picture. 
Eq. (^) can be rewritten as follows 
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y(r,t) 



< $o|t/tft, -oo)[0f + 4>l]U{t, -oo)|$o > 
< $o|C/Hi,-oo)[/(t,-oo)|$o > 



(62) 

where |$o > represents the ground state of an interacting 
electron system described by the free plasmon hamiho- 
nian i/o, and U{ti,t^) is the evolution operator, 



(63) 



T being the chronological operator, and H[/H^^^ repre- 
senting the hamiltonians Hi / H^xt in the interaction pic- 
ture. 



A. Linear Approximation 

Up to first order in the external perturbation one finds, 
after introduction of Eq. ( |63| ) into Eq. (|6^), the linear 
contribution to the induced potential: 

Viir,t)=Vf + 



(64) 



where 



dp 



x^q,p(q- v) e' 



B 

27r ' 

iq-(r||-vt) 



C(^)/q'!p(^0) 



(65) 



and 



/q^W/K^o) 



(2^)2 

xD^(q- v) e"i ("-||-^'*), 



(66) 



D^p{cj) and D^{uj) representing retarded Green's func- 
tions for the operators Xqp{t) and Xq{t), respectively: 



and 



Eq. (p 



2uii 



K.J' 



2-^ 



L0{UJ + i T]) 



(67) 



(68) 



agrees with the linear contribution to the in- 
duced potential obtained, within the semi-classical hy- 
drodynamic model, as the difference between the Fourier 
transform of the total potential of Eq. (|l^) and the ex- 
ternal potential Vcxt{j^,t)- Within the semi-classical ap- 
proach, the role played by bulk and surface plasmons 
goes unnoticed; however, the quantized hydrodynamic 
model provides explicit separate expressions for the con- 
tributions to the induced potential coming from the cou- 
pling with bulk and surface plasmons. The role that bulk 
and surface plasmon excitation plays on the energy loss 
of charged particles interacting with metal surfaces has 
been investigated recently,EEI showing that bulk plasmons 
are excited even in the case of charged particles that do 
not penetrate into the solid. 



B. Quadratic Approximation 



Quantizing the hamiltonian Hi of Eq. (^5|) on the ba- 
sis of both bulk and surface plasmons is a difficult task, 
because of the interaction between bulk and surface plas- 
mon fields. Hence, for the description of the nonlinear 
response to a charged particle moving outside of a semi- 
infinite medium, now we neglect bulk-plasmon contribu- 
tions and find the following expression for the quantized 
Hi hamiltonian: 



H 



l = i^EE[-V,q. Xq,(i)Xq.(i) 

qi q2 

+^q..q. Xi{t)x-IM X^(q,+q,)(i), 



(69) 



where 



7qi+q2 



|qi + q2| 



(7qi - 91 ) (7q2 - 92) 2 "0 91 92 UjS^ ^qi+q2 

|qi +q2|7qi+q2 7qi 7q2 



(gi + 27qi) (92 + 27q2) (|qi + q2| + 27qi+qJ_ 
91 92 (qi • q2 - 7qi7q2) 9l 7q2 (qi ' q2 - 7qi 92) 



7qi + 7q2 + 7qi+q2 7qi + 92 + 7qi+q2 

92 7qi (qi ■ q2 - 91 7q2) 7qi7q2 (qi ' q2 " 91 92) 



91 +7q2 +7qi+q2 



and 

I qi,q2 



1/2 



91 + 92 + 7qi-|-q2 



91 92 |qi +q2|7qi 



(70) 



U)^, UJ^„ LU^ 



7q2 7qi+q2 



qi^q2^qi+q2 + 2 7qi ) 

1/2 



(92 + 27qJ (|qi + qzl 27q,+qJ_ 

(91 + 7qi) (92 + 7q2) (|qi + q2| + 7qi+q2) 
7qi + 7q2 + 7qi+q2 



(71) 



The first term in Eq. ( p^ ) comes from the kinetic en- 
ergy of fluid flow, J dm |V'0p/2, while the second term, 
which is proportional to comes from the internal en- 
ergy G[n\ of Eq. (^). 

Introducing Eqs. (H) and (|9|) into Eq. (|6|), and 
Eq. ( |63| ) into Eq. (62), the quadratic contribution to the 
induced potential is found to be given by the following 
expression: 



d^q d^qi 

(2^(2^2 --q 



V2{r,t)- ^1 

X /|-q, (^o) D^l (^) Dl {ui ) Dl_^^ (a; - ui 



X [UJ UJi A_q^qi -Uj{uj-UJi) A_q^q_qj 

+ a;i(cc>-c^i)Aq„q_q, +3 7'q,_q,] e"!' ("^n , (72) 
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where a; = q • v and wi = qi • v. Linear and quadratic 
contributions to the induced potential (see Eqs. (|6j) and 
([7^)) can be represented diagrammatically as in Fig. 1. 
As for the quadratic contribution, the external perturba- 
tion (white circles) acts twice on the electron gas through 
plasmon propagators (wavy lines), thereby creating an 
induced potential at point r and time t (crosses). 

In particular, for a stationary charged particle we have 
Lu ~ LUi — 0, thereby only the internal kinetic energy of 
Eq. (|^) contributing to the quadratic induced potential. 
For the quadratic contribution to the image potential of 
a stationary charged particle we find 



-Zf 



f3' 



2 , ,2 



6 TT 



dqq dqiqie 



7q 7qi 7q-qi G 



-(q-qi) 



s,,,s 



(7q + 7qi +7q-qi) 



(g + 7q) (qi+7qi) (|q-qi| +7q-qi) 

(9 + 27q) (gi +27qJ (|q-qi| +27| 



(73) 



q-qi; 



which in the case of a non-dispersive electron gas (/? = 0) 
vanishes. Here 9 is the angle between q and qi, and in 
the limit as zq — s- 00 Eq. (|73|) yields 



.3 0-41 

'1 ISttujIz*' 



(74) 



This quadratic contribution to the image potential is half 
the result obtained within the semi-classical hydrody- 
namic model (see Eq. (pO|)). In the case of a stationary 
charged particle the whole quadratic contribution to the 
image potential comes from the second term in Eq.(45), 
i. e., from the linearly induced electron density acting 
twice on the external charge. The total electron density 
ni induced at the surface by a stationary charged parti- 
cle is, in the limit zq — s- 00, times the electron density 
nf induced through coupling with surface plasmons (see 
Appendix A). As a consequence, the total quadratic con- 
tribution to the image potential is in this limit (see Eq. 
(40)) twice as large as the quadratic contribution of Eq. 
(74), which has been deduced by neglecting the coupling 
with bulk plasmons. 



IV. RESULTS 

Fig. 2 shows plots of the linear contribution to the 
image potential of a particle with unit charge e (Zi = 1) 
traveling parallel to the surface of a semi-infinite electron 
gas characterized by a static electron density uq equal to 
the average electron^ensity in the conduction band of 
aluminum {vs ^ 2)xB These plots are shown as a func- 
tion of Zo, the distance from the surface, and the speed 
is taken to be w = (Fig. 2a) and v = 2 (Fig. 2b). Con- 
tributions from couplings with bulk and surface plasmon 



fields, as obtained from Eqs. ( pSj ) and (g^), respectively, 
are represented separately by dashed and dotted lines, 
and the total linear contribution to the image potential, 
obtained from either Eq. ( ^4| ) or Eq. (|6^ ) , is represented 
by a solid line. For comparison, the classical image po- 
tential of Eq. (^ ) is represented by a dashed-dotted line, 
showing that it converges with the full linear result when 
the distance zq is well above the screening length, i. e., 
Zo ^ V^Too for V — and zq ^ 2 y^rjoo for v = 2. 

Quadratic contributions to the image potential of a 
particle with unit charge e {Z^ = 1) traveling, as in Fig. 
2, parallel to the surface of a plane-bounded electron gas 
are depicted in Fig. 3. The electron-density parame- 
ter rg and the velocity of the external charged particle 
are the same as those considered in Fig. 2. In the case 
of a stationary charged particle (Fig. 3a), the contri- 
bution from coupling with the surface plasmon field, as 
obtained from Eq. (O), is represented (dotted line) to- 
gether with the total quadratic contribution (solid line), 
obtained from Eq. (p8|). At large distances from the sur- 
face (zq 00), the total quadratic contribution to the 
image potential is twice as large as the quadratic con- 
tribution from surface plasmon excitation, as discussed 
after Eq. ([7^. At smaller values of zq, the quadratic 
contribution from the bulk channel becomes dominant, 
the total quadratic image potential being near the surface 
larger than the quadratic contribution from the surface 
channel by a factor of ~ 10. The quadratic contribu- 
tion from the surface channel to the image potential of 
a charged particle moving with speed v = 2, as obtained 
from Eq. (|7^), is represented by a dotted line in Fig. 3b. 

Fig. 4 exhibits by a solid line, as a function of the dis- 
tance Zq from the surface, the ratio between full quadratic 
(solid line of Fig. 3a) and linear (solid line of Fig. 2a) 
contributions to the image potential of a stationary parti- 
cle with unit charge e {Zi — 1). For comparison, the ratio 
Zi 1.93 X 10-2 (zo/r^)-^, as obtained from Eq. (|l]), is 
represented by a dashed-dotted line, showing that it con- 
verges with the full ratio (solid line) when the distance 
Zq is well above the screening length, i. e. 



V. CONCLUSIONS 



First of all, we have presented semi-classical and quan- 
tized hydrodynamic models to obtain the quadratic elec- 
tronic response of a semi-infinite electron gas. Then, we 
have derived explicit expressions for the dynamic image 
potential experienced by charged particles traveling par- 
allel to a jellium surface, up to third order in the projec- 
tile charge. In the case of a stationary charged particle 
the total quadratic contribution to the image potential 
has been found to be, at large distances from the sur- 
face, twice as large as the quadratic contribution coming 
from the surface plasmon field. Near the surface, the 
total quadratic contribution to the image potential of a 
charged particle with v = has been found to be larger 
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than the quadratic contribution from the surface channel 
by a factor of ~ 10. As the speed of the moving charged 
particle increases, linear contributions to the image po- 
tential coming from the bulk channel have been found to 
decrease, and quadratic contributions from coupling with 
the bulk plasmon field are also expected to decrease with 
increasing velocity. 

Though nonlinear correctijOas are found to be more im- 
portant far inside the solidc2l than outside, our results 
indicate that the nonlinear image potential is enhanced 
with respect to the linear image potential by a factor that 
is for aluminum as large as ~ 1.15 near the surface in the 
case of a stationary charged particle (w = 0) with unit 
charge e {Zi — 1). At large distances {zq » y/rToo) 
from the surface, the ratio between quadratic and lin- 
ear contributions to the image potential of a station- 
ary charged particle decreases with the distance zq as 
Zi 1.93 X 10~^ {zQ/rs)~^, showing that it vanishes at high 
electron densities. 

As the speed of the moving charged particle increases, 
quadratic contributions to the image potential are found 
to be very small. In particular, in the case of a projectile 
of charge Zi = 10 moving with speed v = 2 near the 
metal surface, contributions to the quadratic image po- 
tential from coupling with the surface plasmon field have 
been found to enhance the linear image potential near 
the surface by a a factor of ~ 1.14. 
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APPENDIX A: 

Here we use both semi-classical and quantized hydro- 
dynamic models to evaluate the first order electron den- 
sity ni, which we assume to be induced by a stationary 
particle (v = 0) of charge Zi that is located far from the 
surface, i. e., zq oo. 

Within the semi-classical hydrodynamic model, this 
quantity is easily found from Eq. (^) to be given by 
the following expression 



(A3) 



the operators Xq,p(^) and Xq(*) being given by Eqs. (|5£ 
and (|5^), respectively, and g^p{z) and g^^z) representing 
bulk and surface coupling functions, 

2^p/(^<,p) 



1/2 



ffq,p(^) 



X {p [ujI + 2 0^ + q^)] cospz 
+qujl sinpz} 



and 



'7 + 7q 



1/2 



(A4) 



(A5) 



The electron density induced by the presence of the 
external perturbing charge is obtained as the expectation 
value of the total electron density operator. Up to first 
order in the external perturbation and in the limit as 
Zq oo, we find 



ni{z) = nf (z) -l-nf(z), 



where 

<(z) - 



,, pWpZ//? 



.2/ 



27I-/3 



and 



nf{z) = Zi 



27r/3 



(A6) 



(A7) 



(A8) 



Eq. (A6) coincides with Eq. (Al), and shows that as 
long as the stationary charged particle is located far from 
the surface the total electron density rii induced at the 
surface (z ~ 0) is, in the limit as zq oo, \/2 times 
the electron density nf induced through coupling of the 
charged particle with the surface plasmon field. 



ni{z) = Zi — 2 

27r/5 Zg 



■e(-z). 



(Al) 



Within the quantized hydrodynamic model, the oper- 
ators corresponding to the induced electron density due 
to bulk/surface plasmons are obtained as follows 



(A2) 



q,p>0 



and 
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FIG. 1. Feynman diagrams representing first (a) and sec- 
ond (b) order contributions to the electric potential induced 
by an external charged particle. The external perturbation is 
represented by white points, and the cross represents a test 
positive unit charge. Wavy lines represent plasmon propaga- 
tors and the black point, joining three plasmon lines, describes 
the nonlinear interaction between three excitations 



FIG. 2. Linear contribution to the image potential of a 
particle with charge Zi = 1 and speed v = (a) and v — 2 (b) 
traveling parallel to the surface of a semi-infinite electron gas 
with rs = 2, as a function of the distance zo from the surface. 
Sofid lines represent the full linear contribution to the image 
potential. Dashed and dotted lines represent contributions 
from the excitation of bulk and surface plasmons, respectively. 
The classical image potential of Eq. (Eq) is represented by a 
dashed-dotted fine. 



FIG. 3. As in Fig. 2, for the quadratic contribution to the 
image potential of a particle with charge Zi = 1. Dotted 
lines represent the quadratic contribution from to the image 
potential coming from surface plasmon excitation. In the case 
of a stationary particle (?; = 0) , the full quadratic contribution 
to the image potential is represented by a solid line. 



FIG. 4. Ratio between quadratic and linear contributions 
to the image potential of a stationary particle {v = 0) 
with charge Zi = 1 that is located outside the surface of 
a semi-infinite electron gas with = 2 at a distance zo 
from the surface. The solid line represents the ratio be- 
tween full calculations of quadratic and linear image poten- 
tials. The dashed-dotted line represents the approximate ratio 
Zi 1.93 X 10"^ {zo/rs)~^, taken from Eq. (^. 
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